Strange nonchaotic attractor in a dynamical system under periodic forcing 
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We observe the occurrence of a strange nonchaotic attractor in a periodically driven two- 
dimensional map, formerly proposed as a neuron model and a sequence generator. We characterize 
this attractor through the study of the Lyapunov exponents, fractal dimension, autocorrelation 
function and power spectrum. The strange nonchaotic attractor in this model is a typical behavior, 
occupying a finite range of the parameter space. 
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I. INTRODUCTION 



The strange nonchaotic attractor (SNA) is an object 
that has chaotic attractor features like fractal dimension 
and nondifferentiability (strangeness) but no exponen- 
tial sensitivity to initial conditions, i.e., its largest Lya- 
punov exponent is nonpositive |Q. It has some analogies 
with trajectories that have been found in a study of the 
Frenkel-Kontorova model || and of the Chirikov- Taylor 
map [^). Aubry has shown Q that the incommensu- 
rate ground-states of the Frenkel-Kontorova model can 
undergo a breaking of analiticity transition, between a 
smooth and a fractal set. Shenker and KadanofT || cal- 
culated the fractal power spectrum of the fractal trajec- 
tory that appears in the Chirikov- Taylor map after the 
breakup of a KAM-like surface a fractal power spectrum. 
This map can be related to an incommensurate driving of 
a nonlinear oscillator. Another example is the accumula- 
tion point of the period-doubling cascade of the logistic 
map Q . However, this attractor occurs in a zero measure 
set in the parameter space. SNA as a typical behavior, 
a finite measure set in the parameter space of a model, 
has been found in the context of nonlinear quasiperiodic 
external forcing (i.e., the forcing of a signal with two in- 
commensurate frequencies) |Q. The study of the SNA 
has also been recently connected with the localization 
problem [g|]. 

A lot of work has been done in order to character- 
ize the features of a SNA; its route of formation, au- 
tocorrelation function and power spectrum. Besides a 
period-doubling cascade, many different routes for the 
formation of an SNA have been proposed: (1) the col- 
lision between a period-doubled torus and its unstable 
parent torus [pTofl ; (2) the progressive fractalization of a 
two-dimensional ergodic torus and, (3) for systems 
with quasiperiodic tori in symmetric invariant subspaces, 
the loss of the transverse stability of a torus [jl2| . Another 
particular feature of the SNA is that its autocorrelation 
function does not decay with the time delay like that of 
the chaotic attrator. It can either be fractal or similar 
to the quasiperiodic case The power spectrum of 

a SNA can be singular continuous and in this case has 



fractal features as discussed in many papers p3|-|18| 

The aim of this paper is to answer the question: are 
SNAs restricted to quasiperiodically driven nonlinear sys- 
tems? We believe that the answer is no. We studied a 
map (hereafter called YOS map) that has been proposed 
in the magnetic context to describe the behavior of an 
analog of the ANNNI model on the Bethe lattice |9]j20| 
and which shows, for external periodic forcing, the oc- 
currence of a SNA pl[ . The same map exhibits typi- 
cal features of a neuron (like activation threshold, nerve 
blocking and rebound behavior) and has been coined as 
a dynamical perceptron of second order (related to the 
dimension two of the map), because it corresponds to a 
two- layer recurrent neural network j2^j. Independently, 
Kanter et al. f23|| proposed a recurrence relation scheme 
known as a sequence generator, which is essentially the 
same map, generalized for any number of dimensions. 
This map can also be viewed as a discrete-time, nonlin- 
ear oscillator, for some values of the parameters. 

Anishchenko et al. have claimed that they found a 
SNA through periodic driving of a map, but Pikovsky et 
al. have shown that it was a chaotic attractor with 
a tiny Lyapunov exponent. We show that Anishchenko 
et aL's attractor is indeed strange chaotic, but ours is 
strange nonchaotic and related to quasiperiodic attrac- 
tors. 

This paper is organized as follows. Section II is ded- 
icated to describing the map and its attractors, mainly 
the strange nonchaotic one. In section III we characterize 
this attractor through its Lyapunov exponents, fractal di- 
mension, autocorrelation function and power spectrum. 
The conclusions are addressed to section IV. 



II. THE MAP AND ITS ATTRACTORS 

The two-dimensional YOS map which this paper is 
concerned with is given by 



x n +i = tanh 

Un+l — X n 



x n - Ky„ + H(n) 



T 



(1) 
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The YOS map was initially proposed |19| to model the 
mean magnetizations (x n , y n ) of the (n th ; (n— 1) )-shells 
of a Bethe lattice, for an analog of the ANNNI model 
p5[-^7f in a constant magnetic field (H was independent 
of n). Here, we extend it for a nonuniform field. T is 
the temperature and n — —J2/J1 where J\(J2) is the ex- 
change coupling between nearest (next-nearest) neighbor 
spins on a Bethe lattice. Yokoi et al. obtained the phase 
diagram of this model at zero field jl9) . Tragtenberg and 
Yokoi studied the effect of finite uniform field P0|] . A 
sinusoidal wave H (n) — H cos(27ru;n) is the particular 
form of H{n) we will adopt throughout this paper, rep- 
resenting a shell-dependent external field/input. 

Kinouchi and Tragtenberg studied the properties 
of the map as a neuron model, where x n is the action po- 
tential of the neuron at time n. 1/T and —k/T are the 
weight factors for the two previous states of the neuron. 
H(n)/T is naturally defined as the external current as a 
function of the discrete time n. They showed that the 
map exhibits many neural features. 

Kanter et al. p3j proposed the following real number 
sequence generator: 



Sj 



tanh 



N 
n=l 



(2) 



and studied this map in the context of time-series and 
neural networks. H(n) could be introduced for repre- 
senting a time dependent input signal. 

From the purely dynamical system point of view, the 
YOS map represents a nonlinear oscillator for the range 
of the parameters considered in this paper, and H(n) 
represents a time-dependent external input. The attrac- 
tors of map YOS can be fixed points, Q-cycles (cycles of 
per iod Q ) , quasiperiodic, chaotic or strange nonchaotic 
For the parameters K = 1,T = 0.5, H = 
the system oscillates with period 6 (see Fig. [ij). But, for 
different values of Hq and oj, keeping the same k and T, 
the attractor becomes richer (see Fig. 0). 



This kind of attractor can also be obtained for other 
values of u>, and is therefore characteristic of this period- 
ically driven nonlinear map. For others examples of SNA 
of this map see [E8| . 
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FIG. 2. First return map of the strange nonchaotic attrac- 
tor of YOS map with k = 1, T = 0.5, H = 0.17 and w = 0.14. 
The initial condition is (xo,yo) = (1, 1). Here are represented 
30 000 iterations, after neglecting the first 10,000. 



III. CHARACTERIZATION OF THE SNA 



A SNA is a fractal object with no exponential sensi- 
tivity to initial conditions (the SNA at the accumulation 
point of the period-doubling bifurcations of the logistic 
map has null Lyapunov exponent but has polynomial sen- 
sitivity to initial conditions p9[ ) . In order to characterize 
this attractor, we investigated the largest Lyapunov ex- 
ponent, fractal dimension, autocorrelation function and 
power spectrum. 
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FIG. 1. First return map for the cycle of period 6 of the 
YOS map with k = 1 , T = 0.5 and H = 0. 



Lyapunov exponents: naive and more accurate 
calculation 



Before calculating the Lyapunov exponent of the at- 
tractor of Fig. ||, let us briefly discuss the sensitivity 
to initial conditions of one of the attractors studied by 
Anishchcnko et al. |18|. 

They proposed a four-dimensional map made up by 
two asymmetrically coupled circle maps, with two cou- 
pling parameters (A and 72). They argued that the for 
72 = the system (1) can be a circle map with quasiperi- 
odic forcing. A small value of 72 will make it an au- 
tonomous four-dimensional map that could show a SNA. 

This four-dimensional map is given by 
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x n+ i = x n + fli - — - sin(27rx„) + 7iy n 
Ztt 

+Acos(27ru„) mod 1, 
Vn+i = liVn ~ tt sin(27rx„), (3) 

ZTT 

K 2 

u n+ i = u n + Vl 2 - — sin(27rii„) + -f 2 (y n + v n ) mod 1, 

ZTT 

K 2 

v n +i = 72 (y n + v n ) - — sin(2nu n ). 

Ztt 



For the set of parameters fii = 0.5, fl 2 = (VZ — 
l)/2,K 2 = 0.03, A = 0.4,7i = 72 = 0.01 and K\ = 
0.8784, Anishchenko et al. claimed the attractor is 
strange nonchaotic. They found a null largest Lyapunov 
exponent within the numerical accuracy of the method 
they used. 

A positive definite largest Lyapunov exponent corre- 
sponds to an exponential expansion of an hypercube of 
initial conditions in at least one of the directions of the 
phase space. In other terms, we can study the sign of 
the largest Lyapunov exponent of a map by studying the 
stretching and contraction of a hypercube of initial con- 
ditions. Here we present a simpler version of this proce- 
dure, studying only the evolution of the distance between 
trajectories generated by only two initial conditions. 

We take two different sets of initial values of (x,y,u,v) 
and calculate the distance d[n) between the trajectories 
generated by each set, as the number of iterations is in- 
creased. That is perhaps the nai'vest way to investigate 
the largest Lyapunov exponent of a map. The first set 
we take as xq = yo — u o = v a = 0.7 and the second as 
x 'a = v'o = u 'a = v a — (0-7 — 10~ 12 ). Figure || represents 
the first 30 000 iterations of the evolution of d(n). 

We can see at first sight that the system is chaotic, 
since the distance between the trajectories with differ- 
ent initial conditions grows exponentially. A simple 
estimation of the slope of the rugged curve leads to 
(0.8 ± 0.3) 10 -3 for the largest Lyapunov exponent A+, 
where we have assumed that the behavior of the distance 
is governed by this exponent and given by 



d(n) ps d(0) exp(A + n). 



(4) 



This result agrees with surprising accuracy with that ob- 
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FIG. 3. Distance between attractors with differ- 

ent initial conditions d(n) as a function of the num- 
ber of iterations n, for the map (|^) for the parameters 

fii = 0.5,n 2 = ^^,K2 = 0.03, A = 0.4,7i = 72 = 0.01 
and Ki — 0.8784. The two set of initial con- 
ditions are xo = yo = ^0 = vq — 0.7 and 
x 'o = Vo — u 'o — v 'o = (0-7 — 10 -12 ). The largest Lyapunov 
exponent is A + = (0.8 ± 0.3) 10" 3 . 



Then, this naive method of checking the sensitivity of 
initial conditions seems to be powerful and we will use it 
as well as the more accurate method due to Eckmann- 
Kamphorst-Ruelle-Ciliberto (EKRC) Q to calculate 
the largest Lyapunov exponent of the attractor of Fig. 



tained by Pikovsky and Feudel_ 
Swinney-Vastano algorithm 



using the Wolf-Swift- 



Fig. |ja exhibits some self-similarity in the behavior of 
d(n) as a function of n. We took many pairs of initial con- 
ditions such that the distances between the initial con- 
ditions from each pair were 10 -2 , 10 -6 , 10~ 10 an<il0 -14 . 
Then we calculated how d(n) vary for each pair as a 
function of the iterations, for the SNA of the Fig. |^. 
The result is represented in Fig. |]b. The various curves 
d(n) x n have a scale invariant like behavior, i.e., for 
various values of the difference in initial conditions the 
evolution with the iterations is rather similar, preserving 
the form for different scales of distance. Moreover, none 
of the curves show exponential divergence, although they 
may vary within few orders of magnitude. It suggests a 
null largest Lyapunov exponent. We confirmed this us- 
ing calculations based on the EKRC method shown be- 
low. This behavior (scale invariance in distance and zero 
largest Lyapunov exponent) is similar to that of a typical 
quasiperiodic attractor. 
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FIG. 4. Distance between trajectories with different initial 
conditions d(n) as a function of the number of iterations n 
for the SNA of the Fig. |^, when: (a) one initial condition is 
xo = Vo = 1 and the other is x' — y' = (1 — 1CP 12 ). No- 
tice the self-similar like structure of the peaks; (b) one initial 
condition is ever taken as xo = yo — 1 and the other has the 
form x ' =y' = (1-A), for A = 10~ 2 , 10" 6 , 10~ 10 and 10~ 14 . 
The distances between the trajectories remain limited even in 
the limit of large number of iterations, pointing out to a null 
largest Lyapunov exponent, and exhibiting a scale invariance. 



Fig. H shows the behavior of the absolute value of the 
largest Lyapunov exponent approximants as a function 
of the number of iterations (neglecting the first 10,000), 
for the attractor of the Fig. ||. This attractor has 
the same shape for many initial conditions: we took 
(x a ,ya) = (±1, ±0.5, 0; ±1, ±0.5, 0). The calculations 
were performed using the Eckmann-Kamphorst-Ruelle- 
Ciliberto method. They do indicate that the largest 
Lyapunov exponent is zero, since the absolute values 
of its approximants scale with the number of iterations 
n as |A+| ~ n -1 . Using the same method, we found 
A_ = -0.709971 ± 0.000001, for the smallest Lyapunov 
exponent. 



FIG. 5. Plot of the absolute values of the largest Lya- 
punov exponent approximants for the SNA of Fig. H as a 
function of the number of iterations n (the first 10 were 
neglected). The initial conditions considered were (2:0,3/0) 
= (±1, ±0.5, 0;±1, ±0.5,0). The absolute value of the ap- 
proximants scales as |A+| ~ n^ 1 , indicating a zero largest 
Lyapunov exponent. The method used here is due to Eck- 
mann-Kamphorst-Ruelle-Ciliberto [31] . 



B. Fractal dimension 

The SNA of the Fig. ^ is a complex geometrical ob- 
ject with a fractal Hausdorff dimension (Dp). In order 
to find out this dimension we have used the box count- 
ing method (3^] . The diagram with the number of boxes 
N(a) visited by the SNA of Fig. ^| as a function of the 
edge length a is shown in Fig.pl The initial condition 
is (xo,yo) — (1, 1), and the first 10 4 iterations were dis- 
carded. We considered the next 10 9 iterations, and box 
edges between a = 10 _1 and a = 10~ 3 . Even with this 
number of iterations, we can see that box edges smaller 
than 10 -2 6 lead to artificially small values of the fractal 
dimension. However, a larger number of iterations were 
computationally prohibitive. The same behavior was ob- 
served in J33[ |. 

Fig. has been constructed by taking ordered sets 
of four consecutive points of Fig. ^| and calculating the 
linear coefficient of the best straight line determined by 
them. Error bars follow from least squares fitting. The 
leftmost point of this figure represents the fitting of the 
four smallest values of log(l/a). The point of order 2 rep- 
resents the fitting of the second to fifth point of Fig. 6, 
counting from the left to the right, and so on. Then, we 
conclude that the fractal dimension is Dp = 1.80 ± 0.09. 
This is the same value found in Ref. |3^] for the attrac- 
tor of Grebogi et al. [Q, but they considered this value 
quite uncertain (since it was obtained from just three 
points). In the same reference, Ding et al. used heuristic 
arguments to conjecture that Dp = 2 and found no con- 
tradiction with the result they numerically found. But 
this value is definitely far from the value we obtained for 



the attractor of the Fig. ||. The evidence points to the 
fractal character of this object. 



J2n=i x(n)x(n + t) 




FIG. 6. Number of visited boxes N(a) as a function of 
the length of the box edge a, for the attractor of the Fig. 



[| The first 
considered. 



10 iterations were neglected and the next 10 



C(r) 



(5) 



The calculation of the autocorrelation function can 
give clues about the nature of the attractor in question. 
Its fractal character can indicate the fractality of the at- 
tractor. However, when we are in the presence of a SNA 
we can observe at least two kinds of autocorrelation func- 
tion: fractal or quasiperiodic |Tq ]. 

Fig. U represents the autocorrelation function of the 
attractor of the Fig. |[ and is very similar to those re- 
lated to quasiperiodic attractors, like that found in Ref. 
|l3f for the strange nonchaotic attractor of the model C 
defined therein. We neglected the first 10 4 iterations and 
considered averages over the next 10 5 iterations. 
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FIG. 8. Autocorrelation as a function of the time delay 
r for the attractor of the Fig. |j| which has the same form 
of the autocorrelation function of a quasiperiodic attractor. 
Here, we discarded the first 10 iterations and considered the 
next 10 5 iterations. 
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FIG. 7. Hausdorff dimension Df of the attractor of Fig. 
^| as a function of the order iof fitting, for sets of four con- 
secutive points of the Fig. |(| The order or fitting i de- 
notes that the slope (and error bars) were calculated for the 
(i, i + 1, i + 2, i + 3) ths points of Fig. ^ (the counting order 
begins from the left to the right). Df is the value of the 
slope of the line determined by these point sets through least 
squares fitting. 



C. Autocorrelation function 



The normalized autocorrelation function of an attrac- 
tor { x n } can be defined as 



D. Power spectrum 

The first step in investigating the power spectrum of an 
attractor given by a sequence {x n } is defining its discrete 
Fourier transform 

N 

s(w, N) = N- 1 ' 2 x n e l2 ™ n . (6) 

ra=l 

Then, we can define the power spectrum of the attractor 
as: 

P(w) = lim < |s(w,7V)| 2 > . (7) 

N— >oc 

The power spectrum of periodic attractors consists of 
(5-peaks at the harmonics of the fundamental frequency, 
whilst in the chaotic case the spectrum is continous. For 
a quasiperiodic case characterized by two incommensu- 
rate frequencies lo\ and t^2, the spectrum contains all the 
frequencies of the form nuj\ + mt02- 
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Many works report p ower spectra of a SNA with singu- 
lar continuous fl5 yL^|l5|]l^ character, like those found in 
some models of quasiperiodic lattices and quasiperiodi- 
cally forced quantum systems [^4|-[3^1 . This spectrum has 
a fractal appearance, where there are peaks weaker than 
(^-functions distributed along a self-similar landscape. 

Fig. ^ shows the power spectrum of the SNA of the 
Fig. y. It has many scales of peaks, exhibiting a fractal 
appearance. We neglected the transient of the first 10 4 
iterations and took the next 10 4 . The detailed study of 
the fractal character of this power spectrum as well as a 
renormalization group approach for it is the subject of a 
forthcoming publication. 
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FIG. 9. Power spectrum of the attractor of Fig. with 
a fractal appearance. The first 10 4 iterations were discarded 
and the following 10 4 considered. Notice the resonance at the 
driving frequency uo — 0.14. 



IV. CONCLUSIONS 

We are able to show that a strange nonchaotic attra- 
tor can result from the dynamics of a periodically driven 
nonlinear oscillator, the YOS map. We have shown that a 
fractal object with zero largest Lyapunov exponent can 
emerge from this dynamics in a finite range of the pa- 
rameters space. The Hausdorff dimension of this object 
is Dp = 1.80±0.09. Its correlation function oscillates like 
those of quasiperiodic attractors and its power spectrum 
has a fractal (or multifractal) appearance. 
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